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Capítulo 3 
El tensor de deformaciones. 


3.1. Deformación en una dimensión. 


Antes de pasar a considerar los casos más generales de deformación en dos y tres dimensiones en el entorno del punto, 
trataremos el caso más sencillo, que realmente ya se ha visto en el tema introductorio, el de la deformación en una dimensión. 


XL Az 
-«l > «a »- 
o nd | | 
| P Q 
10) 
a) Antes de estirar el hilo 
TT =x+s Ax! = Ax + As 
al «al »- 
e . | | 
| | / 
'0) S0 P! Q 


b) Después de estirar el hilo 


Figura 3-1: Deformaciones en una dimensión 


Se considera un hilo elástico que es estirado. Sea O, el origen del sistema de coordenadas, un punto fijo en el espacio (no en 
el hilo, ni fijado al hilo), y P y Q dos puntos fijos sobre el hilo. Al estirar, se produce un desplazamiento s(w) de cada punto, 
pasando los puntos P y Q a las nuevas posiciones P” y Q”, como se muestra en la figura. La nueva posición P” se encuentra a 
x! = 1 +s del origen, la nueva posición Q” se encuentraa x+ Ar +s(2+ Ax) =13+Ar+s(0) + As =3+s+Ar+ As 
del origen y la distancia PQ = Az se transforma en PQ' = Az + As. 


La forma concreta en que cambian estas distancias depende de cómo es el proceso de deformación que se está estudiando. 
Si éste es homogéneo, el comportamiento es del tipo representado en la parte izquierda de la figura 3-2. En este caso, en 
el proceso de estiramiento s es una función lineal de x. La parte derecha de la misma figura muestra un ejemplo de una 
deformación inhomogénea concreta, en la cual s es una función no lineal de x. La forma concreta de la función depende de 
la deformación concreta. 


En realidad, en lo que estamos interesados no es tanto en los desplazamientos absolutos, como en los desplazamientos 


Ss Ss 

$0 So 
O zx O z 
Homogéneo: e uniforme. Inhomogéneo: e no uniforme. 
ses una función lineal de x s es una función no lineal de x 


Figura 3-2: Deformaciones homogéneas e inhomogéneas. En el caso de las deformaciones homogéneas s varía de forma lineal con x, mientras que para las 
deformaciones inhomogéneas esto no es así. 
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relativos, ya que son éstos los que nos darán el estado de deformaciones de un cuerpo. Para ello definimos la deformación 
promedio en el segmento PQ en la forma: 


Alongitud  PQ'-PQ As 


A == 3-1 
longitud original PQ Az en) 


Ex= 


y al considerar la deformación de segmentos cada vez más pequeños, tendremos la deformación en el entorno del punto P 
como: 


, As ds 
o te 6) 


Para el caso de deformación homogénea!, e no depende de x= y por tanto coinciden € y e. En este caso se tiene s = 
sy + ex, donde sy es el desplazamiento respecto del origen del punto que originalmente coincidía con él y e una constante, la 
deformación unitaria, deformación por unidad de longitud medida descontando el efecto de la traslación sin deformación del 
punto P que ya tiene en cuenta sy. La nueva posición de P sería entonces 1' =x+s= sp +(1+e)z 


3.2. Deformación en dos dimensiones. 


3.2.1. El tensor de desplazamientos y el tensor de deformaciones. 


Estudiamos ahora como se produce la deformación de una hoja plana. Como consecuencia del proceso de deformación 
cada punto de la lámina situado inicialmente en una posición 7 = (x,y) sufre un desplazamiento dado por el vector de 
desplazamientos 5. A cada punto le corresponde un vector de desplazamientos y el conjunto de vectores de desplazamientos 
forma un campo vectorial, el campo de los desplazamientos s(r) = $(x,y). Consideramos un punto P de coordenadas 
iniciales 7 y un punto Q infinitamente próximo a P, de coordenadas iniciales 7+ dr (véase la figura 3-3) de forma que dres el 
vector posición de Q medido desde P. Al ser puntos infinitamente próximos nuestros resultados podrán aplicarse directamente 
a deformaciones inhomogéneas. 


Como consecuencia de la deformación el punto P se desplaza a P” (con Pp» = F + 8) y el punto Q se desplaza a Q” (con 
To =TQ+80 = T+dF+35+ds). El vector dr que daba la posición de Q medida desde P pasa a ser d7” = Fq —Fp. = di+d35 
que da la nueva posición relativa de Q” medida desde P”, donde ds da la diferencia entre los vectores desplazamiento de P y 


Q. 
ds = (ds, dsy) 


po Q/ 
| dr = (dx, dy) 
ns diYN | 
Tr =F+8 
di! =dF + d3 Py) 
Q 
ES 
dr 
P(x, y) 
r 
zx 
o 


Figura 3-3: Deformación en dos dimensiones 


Como P y Q están muy próximos y por tanto dx y dy son muy pequeños, podemos obtener la variación de la función de 
dos variables $ = s(x, y) entre dos puntos muy próximos separados dr entre los que la variable x cambia en dx y la variable 


TPara el caso de deformación inhomogénea se tendría la integral s = sy + LL e(a)iduy 2 =xu+s=xu+s0+ Le e(u)dx 
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Figura 3-4: Significado físico de las componentes del tensor de desplazamientos [e;]. 


y en dy 
ds, = 989 dx + PB = Cuad + Cay dy 
Ox Oy dS« Cxa Cay dx ANA" 
== =>  ds=édr, (33) 
9Sy 9Sy dsy Cyx lyy dy 
dsy = oq E + ay Y = Cyadx + eyydy 


donde hemos definido las magnitudes 


a OS» Ml 0Sy _ OSg Ml 0Sy 
NS 5% 


que son cantidades adimensionales pequeñas comparadas con la unidad. 


Como tanto (ds.,, ds) como (dx, dy) son las componentes de dos vectores, un conjunto de componentes [e;¿] que relacio- 
nan las componentes de dos vectores como lo hacen en la ecuación (3-3) deben transformarse como las componentes de un 
tensor de segundo orden. Las e; forman un tensor de segundo orden al que denominamos tensor de desplazamientos. 


Es necesario tener precaución de no incluir en las expresiones de la deformación términos que corresponden a traslación 
como un sólido rígido, ya que éstos dan una variación de la posición absoluta de Q. Como en vez de la posición absoluta de 
Q se está considerando la posición relativa de Q respecto de P (es decir dr) y cuál es la variación de ésta (es decir ds”), se 
está dejando de lado de forma consciente la variación de la posición de P producida no por la deformación entre P y Q, sino 
por la traslación (y rotación) del punto P. 


¿Qué significan físicamente las magnitudes e, que acabamos de definir en (3-4)? Para verlo, estudiaremos qué ocurre 
durante la deformación con dos segmentos dirigidos inicialmente a lo largo de los ejes coordenados. Para el segmento PQ,, de 
longitud dx y dirigido inicialmente según x se cumplirá dy = O mientras que para el segmento PQ,, de longitud dy y dirigido 
inicialmente según y se cumplirá dx = 0 (Véase la figura 3-4). Para PQ,, se tiene 


0S« 
ds. = E dx = exa dx 
(3-5) 
Os 
dsy = y E = Cyadx 


ASÍ, Ca: = ds. /du mide la deformación por unidad de longitud de PQ.,, medida a lo largo de Ox. A su vez, €ya Mide la 
rotación (antihoraria) de PQ.,, ya que el ángulo de giro viene dado por 


2Ón, ey 0 


a => O 


(3-6) 


al ser ds, un diferencial de segundo orden respecto de dx y cumplirse que, para deformaciones pequeñas tan 0 > 0. 


Procediendo del mismo modo para PQ, se obtiene que ey es la deformación por unidad de longitud de PQ, según Oy, y 
que €.» es la rotación en sentido horario de PQ, aP'Q”,. 


Para que este tensor sea una buena medida del estado de deformaciones en el punto P, es evidente que cuando no hay 
deformación todas las e;; deben ser O. Sin embargo, esto no siempre es así. Para ver esto, consideraremos el caso en que el 
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ay 
/ po / 
Qy | Q 
$ 
po Qy dr” d3 
O . 


Figura 3-5: Efecto de la rotación como un sólido rígido 


cuerpo objeto de estudio gira como un sólido rígido un pequeño ángulo $ en sentido antihorario. Cuando se tiene en cuenta 
el significado geométrico de las e;, se tiene un tensor de la forma: 


cosdb—1 —send 
send cosPb—1 


que, en el límite de pequeños ángulos queda 
0 —9 
Cij = (3-8) 
Pb 0 


Así, observamos que aunque no hay una deformación real, las e,¿ no se hacen 0. La razón de esto es que tenemos incluida 
en el tensor el efecto de la rotación como un sólido rígido, y esta contribución no corresponde a una deformación real. Para 
tener un tensor que realmente caracterice el estado de deformaciones en un punto, tendremos que sacar esta contribución. 
Para ello, descomponemos el tensor de segundo rango que hemos obtenido en la suma de dos tensores, uno simétrico y uno 


antisimétrico [e;¿] = [€,¿] + [ws], con las componentes definidas en la forma: 
1 1 
Ej = 5 lei Fez) y wi = 3(eij — ei) (3-9) 


El tensor [e;;], que es simétrico, es el tensor de deformaciones, y es el que da el estado de deformaciones de un cuerpo en un 
punto. 
1 
Exz Exy lxx ¿(xy + yz) 
= (3-10) 


1 
Eyx €yy 3 (€xy + Cyz) Cyy 


Las componentes diagonales son las deformaciones normales al plano correspondiente y dan las deformaciones unitarias 
en las direcciones Ox y Oy, mientras que las componentes no diagonales corresponden a la deformación angular o cortante, 
tangente al plano correspondiente. Dos elementos de longitud paralelos a Ox y Oy antes de la deformación, forman entre sí 
después de ésta un ángulo 5 — 2€,.. Es conveniente recalcar que la deformación por corte €. y es la mitad del cambio de 
ángulo entre los elementos. 


El tensor antisimétrico [w;¿] estará asociado a la rotación del cuerpo como un sólido rígido. De hecho, en la ecuación (3—7), 
el tensor [e;;] que allí aparece es realmente [w;], ya que no hay deformación, y a su vez es [R* — II], donde I es el tensor 
identidad y 7 es la matriz de rotación que da el cambio en las componentes de un vector cuando pasan de representarse en el 
sistema de ejes xyz original a representarse en el sistema de ejes XY Z rotado un ángulo $ respecto del anterior. 


Así, en el caso general se tiene 
ds = [é]dr = [e + úJdr 
dr” =dF+d3= [14 é]dr = [1+€+ú]dr =[R7? + ¿Jar 
que cuando no hay deformación, sino únicamente rotación como sólido rígido, queda 
dr' =R dr > Rdr* = dr 


Se denomina vector deformación al vector £'= [é]u, donde ú es un vector unitario en la dirección de dr. El vector deformación 
así definido está asociado a deformaciones por unidad de longitud en el entorno de un punto, para un segmento infinitesimal 


3.3. Deformación en tres dimensiones. 39 


inicialmente en la dirección dada y por tanto, es una magnitud adimensional. Se denominan componentes intrínsecas del 
vector deformación a sus componentes paralela (€) y perpendicular (e 1) a dr (y por tanto a ú). Nótese que el vector defor- 
mación no es lo mismo que el vector ds (pág. 36) para un vector unitario según dr. En general ds contiene una componente 
de rotación como sólido rígido y además en deformaciones inhomogéneas la contribución de los puntos originalmente según 
uú || dr va cambiando. 


3.2.2. Deformación homogénea en dos dimensiones. 


Durante una deformación homogénea, las e; son constantes a lo largo del cuerpo y se tiene: 


(r;,r5 == x, y) 


si = (so); + Sr eijY; (1,j=x,4) (3-11) 
J 


donde (so); es el desplazamiento del punto original en el eje Or;. Una curva que antes de la deformación es f(x, y) = 0 pasa 
a ser después de la deformación f(x”, y”) = 0, con 
2 =1 +82 = (Soja + U+CxaT + layY 


(3-12) 
y =Y + Sy = (S0)y + Y + CyaT + yyy 


Por ser una sustitución lineal se tiene que, durante una deformación homogénea: 


. Las líneas rectas siguen siendo líneas rectas. 

. Las líneas paralelas siguen manteniéndose paralelas entre sí después de la deformación. 
. Todas las líneas rectas en una dirección se extienden o contraen un mismo factor. 

. Una elipse se transforma en otra elipse, y en particular, un círculo en una elipse. 


Un 


3.3. Deformación en tres dimensiones. 


Para abordar el estudio de las deformaciones en tres dimensiones se sigue un procedimiento completamente similar al 
seguido en dos dimensiones, por lo que no se repetirá. Se definen nueve componentes de un tensor en la forma: 


es — Oi (1,j=x,y,z) 
A Or; Mas Ts = 20,2) 


(3-13) 


donde €, yy, €22 Tepresentan deformaciones unitarias en las direcciones Ox, Oy y Oz, respectivamente, €. y representa la 
rotación alrededor de Oz de un elemento originalmente paralelo a Oy hacia Ox, €,,, representa la rotación alrededor de Oz 
de un elemento paralelo a Ox hacia Oy, y del mismo modo para los demás elementos. Por el mismo procedimiento que en el 
caso bidimensional se obtiene que este tensor no es una buena medida del estado de deformaciones de un cuerpo en un punto, 
que debe estar dado por un tensor simétrico, para evitar la contribución de la rotación del cuerpo como un sólido rígido. Así, 


se define el tensor deformación a partir de [e¿¿] en la forma e; = He; + esi), quedando: 


Cxx >€xy US €yz) Henz alo €zw) 
Exx €xy Exz 
Eya Eyy Eyz | = 5(eyo + Cry) €yy Heya + €24) (3-14) 


Ezx Ezy €zz 


Meza a €xz) >Czy + yz) Czz 


donde las componentes diagonales y no diagonales tienen el mismo significado que en el caso bidimensional, componentes 
normales y cortantes de la deformación. El vector deformación y sus componentes intrínsecas se definen en la misma forma 
que vimos en el caso bidimensional y tienen el mismo significado. 


Como el tensor de deformaciones es un tensor simétrico, existirá siempre un sistema de ejes (ejes principales) en el que 
este tensor es diagonal con autovalores (las deformaciones principales) reales €, €2, €3: 


Exx €xy Exz e 00 
Eya Eyy Eyz | >| 020 (3-15) 
€zx Esy Ezz 0.0 e 


Además de ser los ejes en los que el tensor de tensiones es diagonal, los ejes principales satisfacen la condición de ser tres 
direcciones mutuamente ortogonales que siguen siendo ortogonales entre sí después de la deformación. Nótese que si hay 
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th 


€l 


€3 


€2 


Figura 3-6: Deformación volumétrica. La traza del tensor de deformaciones da la deformación cúbica o variación relativa de volumen. 


componente de rotación como sólido rígido las direcciones principales pueden cambiar a consecuencia de la misma, pero 
siguen siendo mutuamente perpendiculares. 


El tensor de deformaciones debe tener los tres invariantes bajo cambio del sistema de ejes por rotación fija 


1. 


2. 


3: 


3.3.1. 


Traza del tensor o invariante de primer orden, suma de los elementos de la diagonal del tensor, €x.x + €Eyy + €z2 = 
€ +€e + €3. 

La traza del tensor de deformaciones representa la denominada deformación unitaria de volumen o deformación cúbica, 
variación relativa del volumen. Para verlo, consideramos el cubo de arista unidad de la figura 3-6 que se deforma de 
forma homogénea según los ejes principales, con una variación relativa de volumen 


AV 
=V 


donde hemos despreciado términos pequeños de orden superior. La expresión anterior también es cierta para una defor- 
mación inhomogénea tomando un elemento de volumen infinitesimal en el entorno del punto considerado y para cual- 
quier otro sistema de ejes al ser la traza un invariante del tensor. 


0 = (1+e)01+e)(1+e3) -1=e+€.+€3 (3-16) 


Invariante de segundo orden, suma de los determinantes de los menores complementarios de los elementos de la 
diagonal 


Eyy €yz Exx Exz Exx Exy 


Ezy Ezz 


Ezx Ezz Eyzx €yy 


igual a e9€3 + €1€3 + €, €9 en ejes principales. 
Determinante del tensor o invariante de tercer orden 


Exx Exy Exz 
Eyx €yy €yz 
Eza €zy zz 


que en ejes principales es simplemente e; €2€3 


Deformación homogénea en tres dimensiones. 


Se tiene, para este caso de tres dimensiones que, una vez separadas las rotaciones: 


si = (S0)i + Y + e EijYj (3-17) 
J J 


Para separar por completo la parte de sólido rígido (traslación más rotación) se define: 


Si = y EijYPj (3-18) 
J 


Además de las condiciones vistas para la deformación homogénea en dos dimensiones, durante una deformación homogénea 
en tres dimensiones los planos siguen siendo planos y los planos paralelos se mantienen paralelos durante la deformación. 
Asimismo las esferas (un caso particular de elipsoide) se convierten en elipsoides genéricos. 
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3.3.2. Las ecuaciones de compatibilidad. 


Analizaremos en mayor detalle las ecuaciones (3-13) y (3-14). Como tenemos seis e¡; independientes se tendrán seis 
ecuaciones, una para cada componente de la deformación, pero en función de sólo tres componentes del desplazamiento 
Sa, Sy, Sz. Así, SI CONOCEemos S(x, y, z) podremos obtener directamente las deformaciones. 


Un punto de vista alternativo consistiría en considerar las seis deformaciones e;¿ en función de la posición, y a partir de ellas 
intentar obtener 5(x, y, 2). Pero esto nos da un conjunto de seis ecuaciones para el cálculo de las tres incógnitas s,., Sy, 52, que 
en principio sería un sistema de ecuaciones indeterminado ya que tiene más condiciones que variables y por tanto no permitiría 
obtener los desplazamientos a partir del mismo. Sin embargo, esto no ocurriría si existieran relaciones adicionales entre las 
componentes de la deformación que redujeran el número de ecuaciones independientes a tres. Esto significa que en realidad 
las seis componentes independientes de e,; no pueden definirse arbitrariamente en el espacio si queremos encontrar funciones 
únicas y continuas para describir el desplazamiento 3(x, y, z), de forma que el medio se deforme de forma compatible, es 
decir, sin saltos ni solapes. Por este motivo, deben existir al menos tres relaciones adicionales que nos permitan establecer 
una relación biunívoca entre deformaciones y desplazamientos. 


Para obtener estas relaciones comenzaremos con las definiciones de €... Y €yy Vistas al principio de esta sección. Si deriva- 
mos la primera dos veces respecto de y y la segunda dos veces respecto de x, y sumamos los resultados, se tiene 
2 2 3 3 
O Exxz O Eyy  0%8g 09S, 


Oy? * 012  Ox0y? se 0x?0y 


y si, a partir de la definición de €.,y derivamos respecto de x y respecto de y, se tiene 


Perm 1 8 (5 2) 


Dudy ¿dida 0% 


que, comparando con la ecuación anterior, da 


Oe En e 
Tr -2 xy a 
Oy? il da? Ox0y” a 


obteniéndose de forma similar otras dos ecuaciones, 


Oo Pe. ez 
02? i Oy? á 0y0z a 
2 2 2 

19) Exz 0) Exa 99 Exz (3-19c) 


da? 02 — “OrOz 


Partamos ahora de las definiciones de €-:y, Eyz Y Ezu> 


_1/05% , sy). 1/98, 982). 1 [9% , 9 
Su Oy Ox)” te =3NVoz Oy )' == 2 Voz * dx 


Despejando s; y derivando adecuadamente se tiene 


Os Os 0?s, de 0?s 

cd O A 

dy "Ud "Bos 02 0x0% 

OS. ei EN Ps, a 9 Hza Os, 

de “== Oz “Oy  0Ox0y 

Osy 087 _ Os, l ls, D€yz 
> a 7 oa 0 


Si sumamos y derivamos de nuevo respecto de x se tiene 


Dc De 6 Sa 


Ox0z 0x2 dS 0x0y  0x0y0z' 


que comparado con 
ena PS, 


Oydz  Ox0yOz 
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da 
Pena O (O€xy  Otyz  O€za 
Oy0z Ox ( Oz Ox Oy ) : id 
y de forma análoga se obtienen las ecuaciones 
Ve 19) O€xy , Oyz , O6za 
fi EA A -1 
de Oy ( de a ) (3-19) 
e, O (O€xy , O6yz Oz 
d10y Oz ( Oz dx Oy ) (oia: 


que dan un total de seis ecuaciones que deben cumplir las componentes del tensor de deformaciones para que exista una 
solución continua para el campo de desplazamientos. Se les denomina ecuaciones de compatibilidad de Saint-Venant. 


Ahora bien, necesitábamos tres ecuaciones, pero hemos obtenido seis. No debe extrañarnos porque las ecuaciones anteriores 
no son independientes entre sí, y puede verse que son equivalentes a tres ecuaciones independientes de cuarto orden. 


3.3.3. Forma cuadrática asociada a la deformación y elipsoide de deformación. 


El invariante de segundo orden del tensor de deformaciones (forma cuadrática del tensor de deformaciones) tiene la forma 
> EijPiTPj = 1 
0] 


y es idéntica a la forma cuadrática del tensor de desplazamientos 


> eiii = 1 


0) 
como se ve desarrollando su expresión y eliminando términos. Así, ésta última contiene información únicamente sobre la 
parte simétrica. 


Para obtener la deformación unitaria en una dirección arbitraria analizaríamos cómo cambia la componente según esa 
dirección de un segmento unidad Y dirigido originalmente a lo largo de la misma. Al producirse la deformación, tendremos 
en general cambio de longitud y de dirección, de forma que la deformación produce un desplazamiento del extremo en la 


forma 
Asi = > EjjUj 
Jj 
para deformaciones homogéneas, que utilizamos por claridad. Así, la deformación en la dirección original es la componente 


de A3S'en dicha dirección, 
EA = ye As¡u; = SN esj uu; 
1 UN) 
y, en ejes principales, 
EA = E1uí + EQU + €zu3. 
Esto puede verse gráficamente en la figura 3-7, donde la longitud del radio-vector OP en la dirección de UA es 1/y/€ y su 
dirección normal a la cuádrica en el punto P. 


El elipsoide de deformación surge al estudiar como es la deformación de una esfera unidad cuando los calores principales 
del tensor de deformaciones son €, €2, €3: 


A A A A (3-20) 
E (1+e)2 * (1+e2)2 * (1+e3)2 


siendo esta última expresión la del elipsoide de deformación, que no debe confundirse con la cuádrica del tensor de defor- 
maciones, que puede ser un elipsoide real o imaginario o un hiperboloide de una o dos hojas, mientras que como su nombre 
indica, el elipsoide de deformación es siempre un elipsoide. 


En ingeniería se suele utilizar una notación diferente para las deformaciones cortantes, escribiendo e;¿ en la forma: 
1 1 
Exa 3 V10y 3 Y1z 
Ñ € l COn Yij = 264, (321) 
23Vyz €yy 23Vyz Vij = 4€ij 


1 ¿0 
2Vzx 2 V2y Ezz 
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Figura 3-7: Deformación homogénea en tres dimensiones analizada mediante la cuádrica de deformaciones. 


El 


Figura 3-8: Deformación de una esfera en un cortante 


Así, por ejemplo yy, es igual a la variación del ángulo entre dos líneas que estaban dirigidas originalmente según las 
direcciones Ox y Oy. Aunque a menudo se les llama componentes cortantes de la deformación, deformaciones cortantes o 
simplemente cortantes, en sí no son componentes del tensor de deformaciones, ya que a causa del factor 2 no se transforman 
bajo rotación como las componentes de un tensor. Este y es equivalente al ángulo cortante que veíamos en el tema de 
introducción. 


3.3.4. Deformación plana. 


En el caso de la deformación plana, una de las deformaciones principales es O. Un cortante puro es un caso especial de 
deformación plana que tiene la forma, para un cortante en el plano xy: 


0€0 -e00 
e00 |] >45* > 0€0 (3-22) 
000 000 


donde se representa el efecto de rotar el sistema de ejes 45 grados alrededor de z. La traza es O, es decir la deformación cúbica 
es O. Esto se puede apreciar en la figura 3-8, donde se muestra como es la deformación de una esfera en esas condiciones. 
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+ rotación de —y/2 


2 
day € 


Figura 3-9: Deformación cortante vista desde el sistema en el que se saca el tensor de desplazamientos y desde el sistema rotado y/2 en el que se saca el 
tensor de deformaciones. A menudo se denomina cortante simple al primer caso y cortante puro al segundo. El cortante simple es el cortante puro más una 
rotación de y/2. 


Se puede observar como es cero la extensión en las giradas 45”. Nótese también que en general tensión plana y deforma- 
ción plana no son equivalentes, ni siquiera en un material isótropo, un estado de tensión plana puede dar lugar a un estado 
con deformaciones longitudinales en las tres direcciones y en cambio obtenerse un estado de deformación plana mediante 
tensiones longitudinales en las tres direcciones del espacio. 


El caso del cortante simple que se había estudiado en el primer tema corresponde a un deformación cortante pura más una 
rotación de 27, como se puede ver en la figura 3-9. 


3.4. Expansión térmica. 


Asimismo, se puede utilizar el tensor de deformaciones para especificar la deformación resultante del cambio de tempera- 
tura de un material. Cuando se produce un pequeño cambio de temperatura AT, de forma uniforme en todo el cuerpo, aparece 
una deformación homogénea a lo largo del material. En estas condiciones, las componentes de [e5"] son proporcionales a 
AT: 

ej = 04 AT (3-23) 
tér 


¿y ] es un tensor 


siendo las a; los coeficientes de expansión térmica lineal o coeficientes de dilatación lineal. Como [e 
simétrico, [a¿;] también será un tensor simétrico, al ser el anterior multiplicado por una constante. 


tér 


Por esta misma razón, los ejes principales del tensor [a;¿] coincidirán con los del tensor [e;;”]. En ejes principales se tendrá: 


El dérmica) — a AT Elirérmica) — a2AT ESuérmica) — azAT (324) 


siendo los 1, 2, 3 los coeficientes principales de expansión térmica. Por efecto del cambio de temperatura, una esfera del 
material se convierte en un elipsoide con ejes proporcionales a (1 + a1), UU + a) y U + a3). 


La cuádrica para la expansión térmica es: 


Dis Orar =1 en unos ejes generales (825) 


ar? + ar? + Q3rá =1  enlos ejes principales 


El coeficiente de expansión de volumen o coeficiente de dilatación cúbica a, es la traza del tensor [a;¿], que en el sistema 
de ejes principales es A, = (41 + Q2 + Q3 y en un sistema de ejes general 142 €S 0, = Qga + Qyy + 022. Aunque los a; 
suelen ser positivos, hay materiales que los tienen negativos (calcita, berilio, ioduro de plata, etc...). 


En materiales isótropos, que son buena parte de los materiales que interesan en este curso, todas las o,,; son iguales entre 
sí y se tiene +, = 3a,, donde hemos llamado al coeficiente de expansión térmico lineal a, quedando las relaciones entre las 
dimensiones iniciales y finales en la forma: 


l=l(01+ugAD), V=V(+0a,AT) con A, = 3a (3-26) 


como ya se vio anteriormente. 


Capítulo 4 
La ley de Hooke generalizada. 


4.1. El tensor de las constantes elásticas. 


Se ha visto en el primer capítulo la forma que tiene la ley de Hooke para el caso de una tensión uniaxial y la deformación 
unitaria uniaxial asociada, así como para el caso de una tensión cortante y el ángulo cortante de deformación asociado. 


a=Ee 


7=Gy 0 


y es el momento de preguntarse cual es la forma que toma esta ley de Hooke cuando se tiene el efecto combinado de tensiones 
normales y cortantes en el entorno próximo de un punto del material. Hemos visto que para el caso general, el estado de 
tensiones de un cuerpo se describe mediante un tensor de tensiones [7;;], y el estado de deformaciones del mismo mediante 
un tensor [e;¿], y lo que nos estamos preguntando en este momento es cómo y cuál es la relación entre los tensores de tensiones 
y deformaciones. Como los tensores de tensiones y deformaciones son tensores de segundo orden, estarán relacionados por 
un tensor de cuarto orden, que relacionará las componentes de la tensión y de la deformación en la forma: 


Tij = Y Cijurent (42) 
kl 
o a la inversa 
Ej =D Signa Tia (4-3) 
kl 


donde las C;;11 son las constantes elásticas del material y S;;x1 las deformabilidades del material. Al tensor de las C;¿z1 se 
le denomina tensor de las constantes elásticas o tensor elástico. Se denomina a la relación (4-2) (o a la 4-3) ley de Hooke 
generalizada. 


De este modo, se observa que en general, el resultado de aplicar una tensión genérica 7;¡ es una deformación unitaria 
resultante no nula no sólo en la misma componente, sino en principio en todas las demás. Asimismo, en el caso general, una 
componente cualquiera del tensor de deformaciones puede tener su origen en cualquiera de las componentes del tensor de 
tensiones, o en varias de ellas a la vez. Así, por ejemplo, para la componente 9.,;, de la tensión se tiene: 


Cxa = Tax = Crrvater + Coreytry =- Curvztuz + 
+ CraoyaEye + CueyyEyy ml Cuayztyz ale (4-4) 


+ Curozatz E Cureytzy TT Cuarte 


y de forma análoga para las otras ocho componentes del tensor de tensiones. De este modo se observa que la relación entre el 
tensor de tensiones y el de deformaciones es a través de un tensor de cuarto orden y tres dimensiones que en general tendrá 
81 componentes. 


Un análisis de las propiedades de los tensores de tensiones y deformaciones nos muestra que el número de constantes elás- 
ticas independientes es sensiblemente menor ya que muchas de ellas son iguales entre sí, incluso en el caso más desfavorable. 
Para ver esto consideremos en primer lugar la forma que toma la relación entre el tensor de tensiones y el de deformaciones 
para una componente general 7;; del tensor de tensiones: 


Tij = Cijertas + Cijay€xy A Cijyztya pen (4-5) 
Como el tensor de deformaciones es un tensor simétrico, podremos escribir: 
CijayEry ES Cijyatye = (Cijoy + Cijya xy (4-6) 


y análogamente para las demás componentes del mismo tipo. Así, vemos que las componentes Cijx y Cijiz NO aparecen 
nunca de forma independiente, sino siempre como la suma de las dos. De hecho, cuando se miden experimentalmente estas 
constantes elásticas, no se miden cada una por separado, sino siempre se mide la suma de las dos, y lo que afecta al resultado 
no es cada uno de los valores individuales, sino su suma. Así, la relación entre tensiones y deformaciones es la misma para 


45 
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cualesquiera valores individuales de C;¿x1 y Ci¿1 que den el valor correcto de la suma. La elección más sencilla consiste en 
escoger ambas constantes elásticas iguales entre sí, e iguales por tanto cada una de ellas a la semisuma de las dos, a la mitad 
del valor experimental de la suma: 

Cijri = Cit (4—7) 


Consideremos ahora los desarrollos para las dos tensiones 7;¿ y Tj;: 


Tij = ) Cijki€kl 
kl 


(4-8) 
Tji= > Cjikl El 
kl 
Como el tensor de tensiones es simétrico, 7,¿ = T¡¡ y restando las dos expresiones se tiene: 
Tig — Ti =0 = Y “(Cójra — Cin Jena (4=9) 
kl 

que debe ser cero cualquiera que sean las deformaciones unitarias. Para que esto sea así, se debe cumplir: 

Cijri = Cir (4-10) 


Estas dos condiciones deben de cumplirse siempre y reducen el número de constantes elásticas independientes de 81 a 36, 
lo que es lógico, ya que tanto el tensor de tensiones como el de deformaciones tienen sólo seis componentes independientes. 
Antes de tener en cuenta las propiedades de simetría del medio puede verse, mediante consideraciones energéticas, que 
incluso en el caso más desfavorable (la simetría cristalina más baja, sistema triclínico clase 1) el número de constantes 
elásticas independientes es 21 ya que se cumple Ci¿x1 = Chiij. Para simetrías cristalinas más altas el número de constantes 
independientes será menor, tanto menor cuanto mayor es la simetría hasta llegar a las tres constantes independientes que se 
tiene para un cristal cúbico o a las dos que se tienen para un material isótropo. 


4.2. Isotropía. Constantes de Lamé. 


Se dice que un medio es isótropo cuando todas las direcciones son equivalentes en el mismo. En el caso de un material 
isótropo el número de constantes elásticas independientes se reduce drásticamente de 21 a sólo dos. En un medio isótropo 
los ejes principales de los tensores de tensiones y deformaciones coinciden (los ejes de la cuádrica de tensiones coinciden 
con los ejes de la cuádrica de deformaciones). Esto hace que para un cuerpo isótropo sólo haya dos constantes elásticas 
independientes. Por ejemplo, para la tensión normal d.;.,, la ley de Hooke queda: 


Cxxx = Crerater E Caayy[Eyy ES €z2) (4-1 1) 
y para una tensión cortante genérica queda: 
Tij = Coyaytij = (Corre 7 se0y) €ij para 2 Aj (4-12) 


En vez de las constantes elásticas anteriores, en el caso de materiales isótropos se suelen utilizar las constantes de Lamé, 
definidas en la forma: 
A= Crayy 2p = Curve —T zxayy (4-13) 


quedando la relación entre las tensiones y deformaciones normales en la forma: 


Txa = Ména + Eyy + Eyz) + 2U€xa = AO + 2 Eno = [A+ 2) Exa + AEyy + €22) 


Tyy = A€xa + Eyy + €yz) + 2€yy = A0 + 21€ yy = (A + 20)Eyy + Aza + 22) (4-14) 


zz = Ména + Eyy + Eyz) + 21€22 = AO + 2622 = (A+ 2) €z22 + AMéxa + Eyy) 


donde O = €xa + Eyy + €z2 es la deformación cúbica o deformación unitaria de volumen. Para las componentes cortantes se 
tiene una relación de la forma 

Tay = 21€ y Tuz = 2U€xz Tyz = 2f€yz (4-15) 
Ambas expresiones pueden escribirse (con la delta de Dirac 9, = 1 para ¿ = ¡ y 0 en caso contrario y el tensor identidad I 


de componentes 0;¿) como 
Dij= 2pei; + A00s; > 60=2uc+ 01 (4-16) 
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Analizando las relaciones anteriores, se observa que, en el caso isótropo, la relación entre las componentes normales del tensor 
de deformaciones es sólo con las componentes normales del tensor de tensiones, y que la relación entre las componentes 
cortantes del tensor de deformaciones es sólo con la misma componente cortante del tensor de tensiones. Así, para un cuerpo 
isótropo, cada una de las tensiones normales da lugar en principio a deformaciones normales distintas de cero en todas las 
direcciones normales, pero no a deformaciones tangenciales. 


Las constantes de Lamé A y yy así definidas tienen dimensiones de presión como el módulo de Young, las tensiones, o los 
módulos de rigidez o compresibilidad. Para cuerpos muy rígidos, Á y y son prácticamente iguales, mientras que cuanto mas 
plástico es el material menor es u. 


Las constantes elásticas distintas de cero quedan de la forma: 


Cxavmz Cuuvy Ciuzzz 2p +A 
Cueyy Cuyre Cyyzz Cazyy Crezz Cuzoo A (4-17) 
Cyzyz Cozvz Cyayz (Crono NN Caayy) = 2p 


siendo cero todas las otras constantes elásticas que no están relacionadas con estas a través de Ci; = Cijir = Cjini = Cilh- 


4.3. Isotropía. Coeficientes de elasticidad. 


Escribiendo las relaciones anteriores en la forma: 


Txza = (A+ 2M)€xa + AEyy + €22) Fez: =2 toa 


Tyy = (A+ 210)Eyy + Al€za + €22) Tay = 2M€ny (4-18) 


Ozz = (A+ 2M)€22 + Ana + yy) Tyz = 21€ yz 


es posible despejar las deformaciones unitarias. Para ello se utiliza la regla de Cramer quedando para la deformación unitaria 
según x: 


Torio A A 
Tyy (A+2pu) A 
A AR. MAA OR AA OA + 22 4-19 
Ñ ida A z Rar ea an Cu + 022) Pana) 
A A+2p A 
A A A+2pu 
Se definen las cantidades E y Q en la forma: 
pQu +34) 21(2p +34) 
= ——— Q= 4-20 
A+ 2p A ( ) 
aunque en vez de (2 se acostumbra a usar el cociente: 
E A 
== == 4-21 
“0 2(A + y) Ree 


donde E y y son respectivamente el módulo de Young y el coeficiente de Poisson, que ya conocemos. En función de estas 
magnitudes las expresiones quedan, para el caso isótropo en la forma que ya conocemos: 


1 


Exa = E [Ox _ v(Tyy T O zz) 
1 

yy = E [0yy — V(Tzo + 022)] (4-22) 
1 

Ezz = E (022 = V(T za 17 9yy)| 


Al obtener la expresión de la constante de Lamé y se observa que ésta coincide con el módulo de rigidez, quedando las 
relaciones entre las constantes de Lamé y E yu: 
E vE 


E EEE) ne 
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y los elementos distintos de cero del tensor de las constantes elásticas que aparecían en la ecuación (4-17) como: 


E(l-») 
Crrvs = Em E Co HAS 20) 
vE 
Crayy Cyyar Cyyzz Czzyy Cruzz Cizre A (1+DaA—2) (4-24) 
E 
Cusyz Cuzvz Cyayz (Curva > Cxayy) = 2p == ¡e 


Recordemos asimismo que, para el caso de tensión hidrostática, se tenía, para la deformación volumétrica: 


AV g 
V 3e E (4-25) 


donde KC es el módulo de compresibilidad, que tenía la forma: 


0 


K= == (4-26) 


4.4. Energía de deformación. 
4.4.1. Medio isótropo 


Calculemos el trabajo necesario para deformar un elemento de volumen dV = dx dy dz de un material elástico isótropo. 
Como acabamos de ver las tensiones normales están relacionadas únicamente con deformaciones normales y cada una de las 
tensiones cortantes exclusivamente con su deformación cortante asociada. Así, para calcular la energía total de deformación 
calcularemos por separado las contribuciones debidas a deformaciones normales y a deformaciones cortantes. 


Para las deformaciones normales, utilizando la ecuación (4-16), se tiene para las contribuciones a la energía de los trabajos 
de deformación en cada dirección 


D 


v 
1-2 


Wax = Txa dy d2 dex de = 2G e y | AV d€s: 


0 
9Wyy = Tyy dí dz deyy dy = 2G en e | AV deyy 


Ov 
1-2 


0W.. = 02 de dy de;. dz = 2G p + | dV de: 


que sumados dan 


0v 
dWnormales == 2G sde + Eyy AEyy + Ezz dE zz + A) dv, 
donde se ha tenido en cuenta que dd = de... + de€yy + de. al ser O = Ex + Eyy + €z2. Integrando la expresión anterior entre 
que la deformación es nula (de; = 0) y que alcanza su valor final e;¿ se tiene, para el elemento de volumen dV , 


$ 2 2 2 vr á 
Wnormales =G Exa =E Eyy = Exz =+ 1-9y V 
que da una densidad de energía debida a las deformaciones normales, 
e 2 2 2 2 
Unormales = G 2. EN Eyy + Esa + 1= rs | . (4-27) 


Para calcular la contribución de las energías de deformación cortantes, hemos visto que cada una de éstas va por separado, 
por lo que tendremos 

9W.wy = Taydudz dyzydy = Gay dyay AV 

Waz = Tezdady dyz¿dz = Gx 0d Ya AV 

9Wyz = Tyzdxdy dyyzdz = Gyyadyyz aV 
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que, una vez integradas dan, para la contribución a la energía de las deformaciones cortantes 


G 
Weortantes = 5 [en E e E Ajal dv 


siendo la densidad de energía del medio deformado debida a las deformaciones cortantes, 


G 
Ucortantes = 2 [vs + Na =+ Vie] A (4-28) 


La densidad de energía total debida a la deformación queda entonces como 


vV 


—= 2 2 2 
uUu=G Exa + Eyy Y Ea +], 


1 1 1 
04 5Yay Hz Yzz Ty | > (4-29) 


en función de las distintas deformaciones. También se podría haber obtenido la expresión de la densidad de energía en función 
de las tensiones elásticas, partiendo de p.ej, 


1 1 
0 [Tqz — V(0yy + 022)] = E [A + 1)02x — V(Tzx + Oyy + 022)] 
o, en general, 
1 1 
Gi = [A + 1)05 — V(Tza + Oyy + 022) > de¡¡ = E [( + 1v)do — vd(02x + 0yy + 022)] 


que, utilizado en la expresión del cálculo de la energía debida a las deformaciones normales, y sumado a la contribución 
debida a las deformaciones cortantes, daría 


1 v 
U==5 |0og FO yy + Oe TF ss + Oyy + 022) + 27% + 2722 + 2772 | - (4-30) 


donde Oya + yy + 722 €s la traza del tensor de tensiones. 


4.4.2. Caso general 


Para tratar el caso general basta con tener en cuenta que, al estar dentro del rango de la elasticidad lineal, si una fuerza FF es 
de la forma — kx, el trabajo requerido para un desplazamiento x= es kx?/2. Análogamente, la densidad de energía almacenada 
por unidad de volumen en el elemento deformado (igual al trabajo realizado por las fuerzas eternas por unidad de volumen) 
es 

1 
U=>5D Cinacijen (431) 
ijkl 


en función de las componentes del tensor elástico. 


Apéndice A 
Tensores. 


Obviamente, la orientación del sistema de ejes cartesianos en el que describiremos nuestros problemas físicos puede ser 
escogida en forma arbitraria, aunque nos sea más cómodo trabajar en unos sistemas de ejes que en otros. En tanto se trate siste- 
mas de ejes rotados uno respecto del otro con una orientación relativa fija, las propiedades y leyes físicas deben mantenerse 
igual, ya que no pueden depender del sistema de ejes coordenados escogido. 


Para que esto sea así, las componentes que describen el valor de esas magnitudes físicas en un sistema de referencia o en 
otro rotado deben transformarse de una forma determinada y no de otra. De este modo, un tensor estaría formado por un 
conjunto de cantidades que tienen ciertas propiedades de transformación bajo un cambio de base por rotación, resultando 
éstas en una invariancia formal de las leyes físicas que se siguen cumpliendo de la misma forma en la nueva base. Desde 
un punto de vista matemático un tensor sería una aplicación multilineal (lineal en cada uno de sus argumentos), pero nos 
bastará nuestra definición inicial para empezar a trabajar con tensores y ver la relación entre las componentes de la magnitud 
física expresadas en un sistema de referencia ortonormal xyz y las que tendría en otro sistema de referencia ortonormal XYZ 
rotado respecto del anterior. Antes de ello, veremos cómo caracterizar un cambio de sistema de ejes por rotación. 


A.1. Rotación de un sistema de ejes. 


Comenzaremos estudiando la relación entre los vectores unitarios cartesianos en uno y otro sistema de referencia. En 
función de los vectores unitarios del sistema xyz, el vector unitario Ux a lo largo del eje X del sistema rotado se escribe 
como 

UÚx = (Ux Uy )Ug + (Ux Uy) Uy + (UxU¿)Uz 
donde Ux U, representa la proyección del vector unitario dx en la dirección dada por el vector unitario d.,, es decir, la 
componente de ux según U,, Ux Uy SU componente según Uy y UxXU. su componente según 4,. Como todos ellos son 
vectores unitarios, cada una de estas componentes no es más que el coseno del ángulo que forman los dos vectores unitarios, 
que a su vez es uno de los cosenos directores del vector unitario Ux en el sistema de referencia xyz. Análogamente los 
vectores unitarios Uy y Uz se escriben, en función de los vectores unitarios del sistema 1yz, como 


> 


Uy = (Uy Uy) Ug + (Uy Uy) Uy + (Uy U¿)Uz 
uz = (UzÚs) Us + (UzUy)Uy + (úzUz)uz 


Si denominamos ov¿¿ al coseno del ángulo formado por el eje 1 del sistema de ejes rotado y el eje j del sistema de ejes sin 


rotar (a . 1 Y2), podemos resumir todo esto en la tabla (cuyas columnas también valen para obtener ú,, Uy y Uz) 


UX Ox Azxy zz 
Uy | Qyz yy yz 
uz zx Qzy Azz 


que nos da la matriz de rotación 


QOxx Axy Azxz 011 412 413 
=. a XYZ >TY2Z 
R=| Qyx Qyy Aye | = | 21 A 03 con Oj = AJÍ = Us Us (A-1) 
Azx Azy Azz 31 A32 033 


En lo sucesivo utilizaremos indistintamente las notaciones xyz o 123 para denominar los sistemas de ejes y etiquetar las 
componentes de vectores, tensores de segundo orden y matrices de rotación. 


No todos los «¿¿ son independientes entre sí. Por ejemplo, para una rotación en uno de los planos coordenados, el coseno 
director de un eje 2 respecto de un eje j es exactamente menos el del eje ¿ respecto del eje ¿, ya que un ángulo es 180 grados 
menos el otro (véase el pie de la figura A-2), lo que hace que el coseno cambie de signo y la matriz de rotación quede 
antisimétrica en este caso. Sin embargo en el caso más general la matriz de rotación no será ni siquiera antisimétrica, aunque 
sólo tendrá 3 componentes independientes. 


50 
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0(x, Y) = arc cos(Qy) 


0(z, Y) = arccos(Qy2) 0(X, y) = arccoXQ;y) eS 


0(Z,y) = arccos(0zy) 


>= --- > 
- Y >x 
E cosg send 0 
1S R.=| —senó cosó O 
e 0 01 
X Figura A-2: Para una rotación en el plano de la figura, donde se gira un 


ángulo 0 alrededor del eje z en sentido positivo, 0(X, y) es 90 — 6 mien- 

tras que O(x, Y) es 90 + 6. Sumando ambas expresiones se observa que 
Figura A-1: Rotación de un sistema de ejes. Los cosenos de los ángulos  0(x, Y) = 180—0(X, y) y por lo tanto, sus cosenos serán iguales y de signo 
formados por cada uno de los distintos pares de ejes nuevo y antiguo dan los contrario, como se ve en la matriz de rotación para este caso. Sin embargo, 
coeficientes que permiten expresar los vectores unitarios de un sistema de para una rotación genérica en el espacio la matriz de rotación no será antisi- 
coordenadas en función de los vectores unitarios del otro. métrica. 


La matriz de rotación es una matriz unitaria (su determinante es la unidad) de componentes reales y su inversa es su 
transpuesta, 
011 A21 Q31 
RT = Ri = Q12 AU22 QU32 con Oij = UU Í 
013 23 033 


1: sistema rotado 


De Ñ A-2 
j : sistema sin rotar. pal 


La matriz de rotación y sus componentes únicamente dependen de cuáles son los sistemas de ejes que relaciona, y en ningún 
caso del sistema de ejes desde el que estamos observando la transformación. 


A.2. Tensores de orden cero. Escalares. 


Un escalar es un ejemplo de tensor del rango más simple. Cuando una propiedad física está representada por un escalar, el 
valor de éste no depende del sistema de ejes coordenados utilizado. Por ejemplo, la temperatura de un sólido en un punto no 
depende de si tomamos los ejes del sólido de una forma o de si rotamos estos ejes. 


A.3. Tensores de orden uno. Vectores. 


Un tensor de primer orden se representa por un vector. Cuando una magnitud física es de carácter vectorial, su módulo no 
depende de la orientación escogida para los ejes coordenados, pero sí sus componentes según estos ejes. Así, el módulo del 
vector es un invariante frente a la rotación, aunque sus componentes cambien. Si denominamos (%.,, v,, V,) a las componentes 
del vector en el sistema de ejes coordenados xyz y (vx , vy, uz) a las componentes del mismo vector en otro sistema de ejes 
coordenados XY'Z rotado respecto del anterior, en ambos casos tenemos el mismo vector, pero con distintas componentes, 
dependiendo del sistema de ejes en el que lo representamos 


Unyz = VU + Uy Uy Fu¿Us = UxYz = UXUX + Uy Uy + UZUZ. 
Si multiplicamos por el vector unitario Ux se tiene 
ux = (UsUx)Uz + (UyUx)0y + (UsUx )U2 = US cos(1X) + 0 cos(yX) + U¿ cos(2X) = Aga Ug + Any Uy + O92U2 


al ser ortonormales entre sí los vectores unitarios de un sistema de ejes. 
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Un tratamiento análogo puede realizarse para las demás componentes multiplicando por dy y Uz. Escrito en forma matri- 
cial, el resultado queda como: 


cos(1X) cos(yX) cos(2X) 
(vx,vy,uz) = | cos(1Y') cos(yY) cos(2Y) Uy 
cos(1Z) cos(yZ) cos(2Z) Uz 


(A3) 
AXx AXy AXz Vx 
= OY AYy AY z Uy 
AZx AZy AZz Uz 
o escrito en forma más compacta: 
— > > =>) — > ! 
Uxyz =R - Unyz = [045] Unyz O 0'=R-0= 0] 0 O U= ) OipUp (A4) 


donde R es la matriz de rotación cuyas componentes son a, = cos(r¿Or;,) = cos(u 

expresión anterior hemos utilizado la notación con prima para denominar al vector expresado en el sistema de referencia 
XYZ. Ésta también es una notación frecuente, pero debe recordarse que en ambos casos se trata del mismo vector, pero con 
componentes expresadas en distintos sistemas de referencia. 


¿XYZUY2y — ¿XYZ , q ay2 
(au ta) = U -U¿. En la 
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A.4.1. Transformación de un tensor de segundo orden frente a la rotación de ejes. 


Consideremos ahora una propiedad física descrita por un tensor de segundo orden, representado por una matriz, como es el 
caso de la conductividad térmica, eléctrica o hidráulica en el caso general. En estos casos, el tensor (al que denominaremos 
por ejemplo 1) aparece en una expresión de la forma Pr P = q, donde y y son dos vectores relacionados entre sí a través 
del operador T, todos ellos expresados en un sistema de referencia Or. Visto desde un sistema de ejes Or”, relacionado con 
el anterior a través de una matriz de rotación JR, esta expresión queda en la forma T” - p” = q”. Si en la primera expresión 
introducimos la matriz de rotación multiplicada por su inversa, ésta deberá seguir siendo válida: 


T.p=g=T-RT-Rp (AS) 


Si ahora rotamos este vector, se tiene: 
R-¿S3R-T-R*-R-p (A-6) 


Ahora bien, R -G = q' y R-P = P”, por lo que para que ambas expresiones sean equivalentes, cada una en su sistema de 
referencia, deberá cumplirse: Ñ Ñ 
R-T-RI=T' (A-7) 


que es la forma en que se transforman las componentes de un tensor de segundo rango frente a la rotación, 


Es decir, si se tiene un conjunto de nueve cantidades dispuesto con un cierto orden en una matriz 3 x 3, para comprobar si éstas 
son las componentes de un tensor de segundo orden será suficiente comprobar si relacionan dos vectores (cuyas componentes 
se transforman como las componentes de un tensor de primer orden). Esto ocurrirá siempre que los dos vectores representen 
propiedades físicas reales, ya que éstas se transforman de esta forma. 


A.4.2. Ecuación de autovalores. 


Se denominan ejes principales, vectores principales, o autovectores de un tensor al sistema de ejes coordenados en el cual 
el tensor tiene forma diagonal. En el sistema de ejes principales, Or*”, los ejes propios del tensor son obviamente (1, 0, 0), 
(0, 1,0) y (0, 0, 1) y como en este sistema de ejes el tensor es diagonal, claramente se cumple, para el vector propio y,” enel 
sistema de ejes principales ep: 


PE q= (A-8) 
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denominándose a A autovalor, valor propio o valor principal de T con el autovector vu)". Aquí utilizamos el subíndice 
p para denotar que nos referimos al vector principal (o propio) y el superíndice ep para denotar que nos encontramos en 
el sistema de coordenadas principales. Al ser T' un tensor, esta ecuación debe ser cierta en cualquier sistema de ejes. Lo 
podemos comprobar de forma sencilla, considerando otro sistema de ejes que se encuentra rotado respecto del anterior con 
una rotación caracterizada por una matriz de rotación R, si hacemos: 


RTERA RISA RE 
Esto da la misma expresión anterior (A—8) pero en el nuevo sistema de ejes coordenados: 
Td) =A-0, (A-9) 


Así, cuando un vector Y, está dirigido según una dirección principal, el vector resultante de aplicar 7' sobre el vector 1, es 
proporcional al mismo, con independencia de la orientación del sistema de referencia respecto del que se representan T' y 
Up (pero con ambos expresados respecto del mismo sistema). Pasando todo al primer miembro, la ecuación de autovalores 
queda: 

(T — Av, =0 (A-10) 


donde I representa la matriz identidad I,,; = 0,¿. Al ser Y un vector no nulo, debe cumplirse que el determinante de (7 — AI) 
es cero, es decir: 
Ti=A4. Tag T13 


Ta Ta2-A Ta |=0 (A-11) 
T31 Ta  T3z3=A 


A.4.3. Autovalores y autovectores reales para los tensores simétricos de segundo orden. 


Se dice que un tensor de segundo orden es simétrico si sus componentes satisfacen la condición T;, = T;;. La mayor 
parte de los tensores de segundo orden que describen propiedades físicas son simétricos (una excepción es p.ej. el tensor 
termoeléctrico). Un tensor de segundo orden es antisimétrico si sus componentes cumplen la condición T;¿ = —T';i. 


Para un tensor simétrico de segundo orden, existe siempre un sistema de ejes real en el cual el tensor toma la forma 
diagonal, con las tres componentes reales. Si no es simétrico aparecerían valores complejos. Como hemos visto, a los ejes de 
este sistema se les denomina direcciones o ejes principales y a los elementos diagonales del tensor en este sistema de ejes 
valores propios o valores principales del mismo. Cuando alguno de éstos está repetido se dice que hay degeneración. 


Ti Tio Tiz Ti 0 0 T00 
Ta To T33 >R> 0 T3o 0 = 0 Ta 0 (A-12) 
T31 T32 T33 0 0 T3 0 07 


Las propiedades de invariancia bajo rotación del sistema de ejes que obtenemos a continuación las referimos a un tensor 
simétrico de segundo orden, ya que si no, los ejes principales del tensor no serían en general ejes reales. Sin embargo también 
se satisfarían con los valores complejos que se tendrían en el caso de un tensor no simétrico, frente a rotaciones en ejes reales. 


A.4.4. Invariantes bajo rotación de un tensor de segundo orden. 


Para un tensor genérico de segundo orden, la ecuación de autovalores (A—11) queda como: 


AB + (Ti1 + Ta2 + T33)A?— 


(A-13) 
—(Ta2T33 + Ts3Ti1 + Tis Toa — TozT3a — TisT31 — TioTa1) A+ A=0 
que, para el caso de un tensor simétrico, pasa a ser 
A + (Ti + Tao + T33)A*— 
(A-14) 


—(Ta2T33 + T33Ti1 + Ti1 To TES po TEJA | A=0 


de soluciones A = A1,A2,A3 = 71, T», T3, todas ellas reales. En la ecuación anterior se ha utilizado A para el determinante 
de [T;,|, que para un tensor genérico de segundo orden como [T;¿] tiene la forma general. 


A = Ti Tr T33 + Ti2Ti3T31 + Ta1T32 713 — TigTo2T31 — Tio To1T33 — TsoTo3Ti1 (A-15) 
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Si el sistema de ejes está dirigido a lo largo de las direcciones principales, los elementos de fuera de la diagonal son todos 
nulos, y el determinante queda A = A14243 = 7, T73T3, donde A1,A2,A3(= T1, Ta, T3) son los valores propios del tensor 
obtenidos a partir de la ecuación de autovalores (A—11). 


Las soluciones A; de la ecuación (A—14) no dependen del sistema de referencia, y por lo tanto, los coeficientes de A en dicha 
ecuación tampoco pueden depender del sistema de referencia. Se dice que son invariantes frente a la rotación. Se denomina 
traza a la suma de los elementos de la diagonal. En un tensor simétrico de segundo orden, la traza o invariante de primer 
orden es un invariante frente a la rotación, es decir: 


L = NT = Y Ti (A-16) 


Existen además otros dos invariantes para un tensor simétrico de segundo orden, el denominado invariante de segundo orden, 


—=Lo = —(To2Ts3 + Ts Ti + Ti Tos — Tag T3o — TigT31 — TioTo1) 


(A-17) 
= (TT, + TaT3 + T3T1) 
que también puede escribirse como la suma de los menores complementarios de los elementos de la diagonal, 
T22 Tas Tas Tis lu Ti 
L = A-18 
a | T32 T33 Ta Ti Ta Toa ( ) 
y el determinante o invariante de tercer orden, 

lí =A=T BT; (A-19) 


correspondientes a los coeficientes de los términos lineal en A e independiente. Todos ellos permanecen invariantes bajo 
rotación del sistema de ejes. 


A.4.5. Rotación de ejes en un plano principal: el círculo de Mohr. 
El círculo de Mohr constituye una representación gráfica muy sencilla y útil para estudiar el cambio de las componentes de 
un tensor de segundo orden cuando se produce una rotación del sistema de ejes en uno de los planos principales. 


Consideremos de partida un determinado tensor de segundo orden en su correspondiente sistema de ejes principales, con 
valores principales 01, 92 y 93, ordenados de menor a mayor. 


Ñ 01 0 0 
TP?=| 002 0 (A-20) 
0.0 93 


Si se produce una rotación de ejes en uno de los planos principales, por ejemplo, el plano xy (rotación alrededor del eje z) 
caracterizada por una matriz de rotación 
cosg senó O 
R= | —senó9 cos9 0 (A21) 
0 o 1 


las componentes del tensor se transformarán bajo la rotación como RT'R”!, es decir: 


cosó sen0 0 0100 cosO —senó 0 
—senó coso O | - 0097 0 |-| senó coso 0 
0 O 1 0 003 0 0 1 
que desarrollado queda 
01 cos?0 + 07 sen?0 —o¡sen0cosó + oasenOcosó O 
T"=RT*PR= | —0,sen0 cosó + 0, sen Ó cos O 0 sen?0 + 07,c08?*0 0 (A-22) 
0 0 03 


Es decir las componentes diagonales 911 y 922 del tensor en el nuevo sistema de ejes son 


011 =01,c08?0 + 02 sen?8 (A-—23a) 
077 =0,sen?0 + 02c08?8 (A-23b) 
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y las no diagonales 7,2 (iguales a 721 para un tensor simétrico) quedan 


T129 =—01sen0cosÓ + 02 sen cosO (A-23c) 


Demos a estas expresiones unos pequeños retoques. Por ejemplo a la (A-23a), 
— 2 29 — 
011 = 01 COS 0 + 02 sen 0 = 


(cos? + cos*8) + S(sen?9 + sen?0) = 


(cos? +1-sen*0) + 2 (sen?9 +1 -— cos?9) = 
SE (cos? —sen*0) — (cos? — sen 0) 


2 2 


Con retoques similares en la ecuación (A-23b) y recordando que sen 29 = 2 sen 0 cos 0 y cos 20 = cos? 9 — sen? Y se pueden 
reescribir las ecuaciones anteriores como 


Pd - de o a A cos 29 (A24a) 
is - $ => cos 29 (A-24b) 
ma = 22 > 2 sen 20 (A-24c) 


01 + 02 
2 


01 


e Pa ] . Ln DO 
Pero esta es la ecuación paramétrica de una circunferencia centrada en el punto ( 5 o) y de radio A como 


la que se muestra en la figura A-3 y a la que se denomina círculo de Mohr. 


Esto mismo se podría haber hecho para cualquiera de los otros dos planos principales, teniéndose distintos círculos de Mohr 
para cada plano principal como se muestra en la figura A-4, 


El círculo de Mohr permite obtener de forma muy sencilla las componentes de un tensor cuando se produce una rotación 
de un ángulo 0 en un plano principal. Para ello, si el ángulo realmente girado es 0, giramos 20 el segmento AB en torno al 
origen del círculo de Mohr. Los puntos de corte de este segmento con el círculo de Mohr dan las componentes del tensor para 
este giro. Así, las proyecciones de estos puntos sobre el eje horizontal dan las componentes diagonales transformadas, 011 y 
022 mientras que la proyección sobre el eje vertical del punto B da la componente no diagonal transformada, de acuerdo con 
las ecuaciones (A—24). En cada uno de los planos principales el máximo valor de la componente no diagonal tendrá lugar 
cuando AB esté vertical, es decir, cuando 0 = 45% que dará 29 = 90”. Obviamente, como se ve en la figura A4 el máximo 
valor no diagonal tendrá lugar en el plano en el cual los valores principales correspondientes estén más separados entre sí, en 


el caso de la figura, el plano [13] a 45” de los ejes principales girando en torno al otro eje principal. 


El círculo de Mohr también permite, una vez conocido que el sistema está en uno de sus planos principales, obtener valores 
y direcciones principales a partir de cálculos trigonométricos sencillos 


A.4.6. La superficie indicadora o forma cuadrática asociada a un tensor simétrico de segundo orden 


Una forma cuadrática es una generalización del producto escalar de dos vectores p' y Ya través de un tensor de segundo 
orden T;;, 


pTI=) pi Taj =1 
7 


Consideremos una matriz 7; de nueve componentes a la que aplicamos dos vectores posición por la izquierda y por la 
derecha de forma que el resultado sea la unidad, 


) a =1= (r],12,13) [011 %2 013 ri 
13 Q91 022 023 ra 
031 Aza U33 P3 


que desarrollado queda 


Duri+Tisrira + Tigrira + Tarrari + Taor3 + Togrora + T3irá + T3argra + Tagrá = 1. 
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Figura A-3: Círculo de Mohr en un plano principal. Nótese que la rotación en el círculo de Mohr, 20 es el doble de la rotación 0 en el sistema de ejes real. 
Los sentidos de ambas rotaciones son idénticos, positivo como se muestra en la figura. 
T 


La 


Figura A-4: Círculo de Mohr en los tres planos principales. El máximo valor posible de las componente no diagonales se produce en el plano principal 
correspondiente a los autovalores extremos. 


Si consideramos el caso de una matriz simétrica, T,¿ = T;¡, esta suma puede escribirse como 
2 2 2 
Ti1rí + Ta273 + T33r3 + 2T13r173 + 2T33r9713 + 2T33r913 = 1, 
que es la ecuación de una superficie de segundo orden (una cuádrica) centrada en el origen y que puede ser bien un elipsoide 
o bien un hiperboloide, según los signos de 7, Ta y T3. 


Si ahora transformamos las componentes de los vectores 7” mediante una rotación del sistema de ejes dada por una matriz 
de rotación R = [a4;], se tiene 1; = Y y AxiTy Y 1 = Y» 037] (nótese que es la transformación inversa), y la ecuación de 
la forma cuadrática queda como 


e / cn / 
> Tij¡Q4i0151,1, =1= ) Tr, ¡P Y] con Te = > Oki 015 Ti; 
ijkl kl ij 


A.4. Tensores de segundo orden. Representación matricial. 57 


(a) (b) (c) 


Figura A-S: Distintos tipos de cuádricas. Cuando 71, T2 y T3 son positivos se tiene el elipsoide de (a), cuando dos son positivos y uno negativo (en el caso 
concreto de la figura, el negativo es 73) se obtiene un hiperboloide de una hoja como (b) y si dos son negativos y uno positivo (en el caso de la figura, el 
positivo es 73) un hiperboloide de dos hojas como en (c). 


que coincide con la forma en que se transforman las componentes de un tensor simétrico de segundo orden. De esta forma 
vemos que para ver cómo se transforman las componentes de un tensor bajo rotación del sistema de ejes basta con ver como 
se transforma la correspondiente cuádrica. 


Una superficie cuádrica de esta forma tiene ejes principales, es decir, unos ejes en los que la cuádrica se representa como 
2 2 2 
Tirí + Toarz + T3r3 = il 


donde los 7, son los valores principales del tensor correspondiente. Si comparamos con la ecuación general de una cuádrica, 


q? y? 2? 


2 ara 


vemos que los semiejes correspondientes son precisamente 1/yT1, 1/WT3 y 1/VT3. 


Si los tres valores 71, T2 y T3 son positivos la superficie es un elipsoide real, y si los tres son negativos tendremos un 
elipsoide imaginario. Si dos de los valores son positivos y uno negativo la cuádrica es un hiperboloide de una hoja y si uno es 
positivo y dos negativos es un hiperboloide de dos hojas como se ve en la figura A-5. 


Veamos una propiedad interesante basada en las características de una elipse, 


Si T,¿ son las componentes de un tensor simétrico de segundo orden que relaciona los vectores Py en la forma 


Pp = T( (y por tanto, p; = T;¡qs) la dirección de p para un q dado puede obtenerse dibujando el vector 
posición 7 de un punto Q de la cuádrica tal que 7'sea paralelo a q y encontrando la normal a la cuádrica en Q. 


Para demostrar esto partimos del tensor T' en ejes principales. El vector %¿ tiene componentes (u;, us, uz) y, si OQ tiene 
longitud r, el vector posición de Q, F tendrá como coordenadas r(u;, uz, uz) = (11,72, r3). Como Q está en la superficie sus 
coordenadas cumplirán T,r7 + Tará + Tari. 


En una cuádrica, si los cosenos directores del vector posición de un punto de la misma desde el origen son (u;, uz, uz) el 
vector (u¡T|, uzT», uzT3) es normal a la superficie de la cuádrica en ese punto. Como 


p= [T.jJd = q(u Ti, uz T>, ugT3) 
2 


A 


Figura A-6: Forma cuadrática asociada al tensor de tensiones. 
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Imaginario 


[T,] 2? — [Tal y? =-—1 


Figura A—7: Trabajando con un hiperboloide real e imaginario en el plano 12 cuando 7, > 0 y T2 < O. El hiperboloide correspondiente es el marcado 
en negro y cuando la dirección deseada lo corta se procede como de costumbre. Cuando la dirección deseada no pasa por el hiperboloide real es necesario 
trabajar con el hiperboloide imaginario complementario (en gris y a trazos), de forma que se busca el punto de corte, se traza la normal y se le da a ésta el 
sentido contrario que en el caso del hiperboloide real. 


p deberá será normal a la superficie en ese punto. 


Esto se ve muy claro en el caso de un elipsoide real, pero si tenemos los tres valores negativos lo que tendremos en realidad 
es un elipsoide imaginario. En la práctica esto no añade gran complicación, se trabaja con el elipsoide real que correspondería 
al (a) cambiado de signo (Tx? + Tay? + T3z2? = —1, que sería un elipsoide real al ser T,,T3,T3 < 0) y se cambiaría de 
signo el resultado, que apuntaría en el sentido contrario, hacia dentro del elipsoide. 


Asimismo, cuando tenemos hiperboloides, de una o dos caras, hay direcciones en las que el radio vector no corta el 
hiperboloide real, por ejemplo, si en la figura A—5-(b) la dirección es próxima al eje 3, o si en la figura A—5-(c) es próxima al 
plano 12. En estas direcciones, el punto de corte sería imaginario, pero para trabajar con él utilizaríamos el hiperboloide real 
complementario (el que se obtendría cambiando de signo al 1, por ejemplo en la figura AS (b) y (c) serían complementarios) 
y trabajaríamos con él en esta zona, invirtiendo el sentido del resultado. Para mayor claridad esto se representa en la figura A— 
7, para el caso de un tensor plano con 71 > 0 y T2 < 0. 


A.5. Tensores de orden superior. 


Además de los tensores que hemos visto hasta ahora, pueden aparecer tensores de orden superior en las leyes de la Física. 
Por ejemplo, el tensor piezoeléctrico es un tensor de orden 3 que relaciona un vector (la polarización) con un tensor (el tensor 
de tensiones) en la forma 


[P,] = [d;51] [05]. 
Las 27 componentes de éste tensor de orden tres se transformarían, bajo rotación del sistema de ejes, en la forma 
di ¡pQj d 
ik = A ipQjq AkrUpgr 
pqr 
la relación entre el tensor de tensiones y el de deformaciones viene dada a través de un tensor de orden 4, el tensor de las 
constantes elásticas, 
[0,5] = [Cigra]l€n] 


que tendría 8l componentes (afortunadamente ni mucho menos son todas independientes) que se transformarían en la forma 


lA 
¿gl = o VipAjqUkr Qs Cpqrs- 
pqrs 


Apéndice B 


La función de Airy para un estado de tensión 
plana. 


Resolver un problema elástico con condiciones de contorno es muy difícil en el caso general, y de hecho durante muchos 
años, algunos grandes matemáticos han trabajado en el análisis de los problemas con condiciones de contorno. Se han desa- 
rrollado una serie de técnicas, algunas muy elegantes, y aunque ninguna de ellas es completamente general, han permitido 
obtener una información preciosa sobre el comportamiento de los materiales deformables. En la actualidad, este tipo de 
problemas se suelen abordar de forma numérica mediante cálculos de elementos finitos, pero algunas de las soluciones exac- 
tas y los métodos para obtenerlas siguen siendo importantes. 


La dificultad de un cálculo exacto estriba en que, cuando se quiere resolver un problema general en elasticidad lineal, 
incluso en el caso más simple de medios isótropos, es necesario tener en cuenta las ecuaciones de equilibrio (2-3), las 
definiciones de las componentes del tensor de deformaciones en función del desplazamiento (3-14) y las relaciones tensión 
deformación para un medio isótropo dadas por (4-22), sin contar con la posibilidad de que simultáneamente se produzca 
expansión térmica, todo ello basado en el conocimiento que tenemos de lo que ocurre en los límites del material tensionado. 
Esto implica ecuaciones diferenciales en derivadas parciales en las que intervienen vectores y tensores, lo que en la práctica 
multiplica el número de ecuaciones. 


Una posible aproximación para intentar simplificar el número de ecuaciones involucradas consistiría en intentar reducirlas 
a una única ecuación diferencial en derivadas parciales para una función escalar, de la que se derivarían las tensiones y 
deformaciones. Esto no puede hacerse para el caso general de tres dimensiones, pero la función de Airy permitirá hacerlo en 
el caso de un estado de tensión plana. 


Supongamos un estado de tensión plano, en el que p.ej., las caras z están libres de tensiones (7, = 0, para 2 = x, y, 2). En 
estas condiciones las ecuaciones de equilibrio (2-3) que deben cumplir las tensiones elásticas se quedan en 


Osa OTxy 


q —= 0 
Ox + Oy +P9 
OTxy 00 yy 

=0. 


Como una forma de simplificar la búsqueda de funciones cuyas derivadas satisfagan las ecuaciones anteriores, G.B. Airy 
definió una función escalar de tensión p(x, y) a partir de la cual se pueden obtener las componentes del tensor de tensiones 
en la forma 


odo) 
Cxrx = Oy? + Q 
9? 
0 =33+2 (B-1) 
__Y 
e OyOx 


que, como puede comprobarse por sustitución, satisface las ecuaciones de equilibrio si la fuerza de masa admite un potencial? 
de forma que (2 cumpla 


00 
x > — —_— > — — B-2 
pg da P9y y (B-2) 


donde 9. y 9y representan componentes de la aceleración debida a fuerzas de masa y por tanto, P9x y p9y son las compo- 
nentes de la fuerza por unidad de volumen. Se denomina a la función ¿(x, y) función de Airy o función de tensión de Airy. 
Nótese que si la fuerza de masa admite un potencial, pero no se satisface (B-2) no es posible definir una función de Airy. 
Afortunadamente, para el caso de la fuerza gravitatoria la ecuación (B-2) se satisface con (9, 9y) las componentes de la 
aceleración gravitatoria en el sistema de ejes considerado, quedando la densidad de energía potencial Q en la forma 


(2 = —PZYa — PY9y- (B-3) 


T Aquí (2 sería realmente una energía potencial por unidad de volumen, una densidad de energía potencial. 
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Consideremos ahora la condición de compatibilidad (3-19a) y las relaciones (4-22) entre tensiones y deformaciones para 
un medio isótropo. Sustituyendo los valores de las deformaciones dados por (4-22) en (3-19a) y particularizando al caso de 
estado plano de tensiones, se tiene 


e 9 Pr 
dy? (Ozx — VOyy) + Da? (Tyy — VO zx) =2(1 +1) 920y 


de forma que las tensiones que satisfacen esta relación garantizan la continuidad de la deformación. Si sustituimos en esta 
ecuación los valores dados en (B—1) se tiene 


op 0% 050) 


a A AT E A 
+ aa + da20y? 


ES 2 2 _ 4 
dt * By =V"(V*p) = VU. (B-4) 


Así, toda función Q(x, y) que satisfaga esta ecuación cumplirá automáticamente las condiciones de equilibrio, de compatibili- 
dad y las relaciones de la elasticidad lineal para un medio isótropo. Hay muchas funciones que satisfacen la ecuación anterior, 
pero para que sea solución del problema objeto de estudio es necesario además que satisfaga las condiciones de contorno de 
tensión y deformación en los límites del medio deformado, lo que complica mucho las cosas. Puede demostrarse que si se 
encuentra una solución que satisfaga simultáneamente las ecuaciones de compatibilidad y las condiciones de contorno, ésta 
es única y por tanto es la solución correcta. 


La función de Airy puede utilizarse únicamente en el caso de tensión bidimensional y cuando la fuerza de masa admite un 
potencial en la forma mostrada en la ecuación (B-2). 


B.1. Función de Airy para un voladizo de sección rectangular sometido a su propio peso 


Consideremos el caso de un voladizo de sección rectangular para ofrecer un ejemplo de obtención de la función de Airy. 
Partiremos de un desarrollo polinómico de quinto orden en x e y, tomadas como se indica en la figura B-1, 


p(x,y) = Asor? + Agz*y + Agoaóy? + Azgga?y? + Ayazy* + Aosyó+ 


+4Agor* + Azraiy + Az22?y? + Arzuy? + Avay*+ 
+Agot* + Agia*y + Ajay? + Aosy*+ 
+A2o2? + Ar1ty + Apay? 


(B-5) 


y tendremos en cuanta que las aceleraciones debidas a fuerzas de masa son 9, = 0 y gy = 9 (nótese que el eje y apunta hacia 
abajo). Así, a partir del desarrollo anterior, las expresiones (B—1) nos dan para Oy:.x, Oyy Y Tuy 


9? 9? 

Oo = ya — P29w — PYY = ya — P9Y= (B-6a) 
24391? + 64o32*y + 1241404? + 20A05y*+ 

+4g21? + 6A13xy + 12404y*+ 

+24123 + 6A03y + 2402 — pgy 


y? 06 

Pq2 7 PE9e — PUYy = az — PIU = (B-6b) 
20Asp1? + 124412*%y + 6Az21y? + 2A93y*+ 

+124A4p2? + 643,1y + 2A99y2+ 

+6A307 + 2491Y + 2490 —= Pg9Y 


E 


y 


Figura B-1: Obtención de la función de Airy para una ménsula sometida a su propio peso. 
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Tay = == 4A412% + 6A3o1*y + 6Ao312y? + 4A1ayó+ (B-6c) 
Ox0y 


+3A311? + 4A921Yy + 3A13y2+ 
+24911 + 24194 + Ar1 


Condiciones en las tensiones cortantes. Todas las superficies externas están libres de fuerzas cortantes. Por tanto, las 
tensiones cortantes deberán ser nulas en todos los puntos de las mismas. 


Para la cara extrema a x = L debe cumplirse 7,y(L, y) = 0 y por tanto 
Tay(L, y) =0= 44144? + (64231 + 3413)y?+ 
+(64321? + 4A09L + 2412)y+ 
+44411% + 34311? + 2491L+A11 


también deberá cumplirse para todo valor de y, por lo que los coeficientes de las distintas potencias de y deben ser todos 0, 


Aja=0 (B-7a) 

2423L + A13 =0 (B-7b) 

34321? + 24221 + Aj2 =0 (B-7c) 

4A4 1% +34311? + 24911 + A11 =0 (B-7d) 


Para las caras superior e inferior, las tensiones cortantes también serán nulas en todos los puntos, 7; (2, +a) = 0. Así, 
utilizando la ecuación (B-6c) se cumplirá 


Tay (1, +a) =0= 4A412* + (464320 + 3431)0%+ 
(6 Az3a? pa 4A924 + 2491)1+ 
+4 440? + 3A13a? t 24124 + Aj1) 


OS 


para todo valor de x, donde el signo superior se refiere a la cara inferior y el inferior a la superior. Todos los coeficientes de 
las distintas potencias de x deben ser 0, y esto debe ser así tanto para la cara inferior como para la superior. Si tomamos, por 
ejemplo, el coeficiente de 2? para las caras superior e inferior y sumamos o restamos ambos el resultado seguirá siendo nulo, 
pero además observamos que la suma de las partes que no cambian de signo debe ser nula, lo mismo que la suma de las partes 
que cambian de signo Sen el fondo es lo que obtenemos al sumar o restar los coeficientes de las caras superior e inferior). 
Así, obtenemos, a partir de la condición de contorno de las tensiones cortantes en las caras superior e inferior, 


As =0 (B-8a) 

As2=0 (B-8b) 

Asi=0 (B-8c) 

Ag =0 (B-8d) 

3 A934? + Ao; =0 (B-8e) 
4A14a7 +24190=0 (B-8f) 
3A130? + Ay1 =0. (B-8g) 


Ahora bien, vimos en (B-7a) que A14 es cero, por lo que, a partir de (B-8e), A12 también debe serlo. Resumiendo y simpli- 
ficando los resultados obtenidos a partir de las condiciones de contorno para las tensiones cortantes, 


A12 = Aja = Az = Az1 = Aza = Agr =0 (B-9a) 
2423L + A13 =0 (B-9b) 

24211 + 411 =0 (B-9c) 

34230? + Ag =0 (B-9d) 

3A13a? + Aj1 =0. (B-9e) 


Condiciones en las tensiones normales. Veamos ahora qué información tenemos acerca de las tensiones normales en los 
límites. 
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Por una parte, la tensión normal de tracción en la dirección x debe ser nula en todos los puntos del extremo libre, 
OxalL, y) = 20Ao0sy* + 124044? + (64231? + 6413L + 6Ao3 == p9)y + 2402 


donde ya hemos tenido en cuenta que 414 = A22 = Aza = 0. Como esta expresión debe cumplirse para cualquier valor de 
y, los coeficientes de sus potencias serán todos nulos, quedando 


As =Ú (B-10a) 
Ap=0 (B-10b) 
64231? + 6413L+6403 —pg=0 (B-10c) 
Aga =0. (B-10d) 


Veamos otra condición más. Tanto en la cara superior como en la inferior, la tensión normal 0, (uv, +a) debe ser nula, lo que 
da las ecuaciones 


Oyyl2, +a) = 20A501* + 124401? + 6Azyz + (+2 4234? E 249,4 + 240 + pga) 


donde se ha tenido en cuenta que 422 = A31 = A32 = A41 = 0 Los coeficientes de las potencias de x deben ser nulos y, 
dentro de éstos la suma de las partes que cambian de signo y la suma de las partes que no cambian de signo, dando 


As =10 (B-11a) 
Ago =0 (B-11b) 
Azo =0 (B-110c) 
As =0 (B-11d) 
24234? + 24910 — pga=0 (B-11e) 


Combinando todos los elementos. Una vez aplicadas las condiciones de contorno anteriores nos queda una sistema de seis 
ecuaciones 


2A93L + A13=0 
2421L+411=0 

34930? + Ag, =0 

341307 +41, =0 

64231? + 6413L + 6403 — pg=0 
249347 + 24910 — pga=0 


con cinco incógnitas. En realidad, las cuatro primeras ecuaciones no son linealmente independientes, luego lo que tenemos 
es un sistema de cinco ecuaciones con las cinco incógnitas que debemos extraer. Esto es muy sencillo, ya que entre las 
ecuaciones (B-9b) y (B-1 le) es trivial despejar A21 y A23 y sacar entonces 411 de (B-9c) y A13 de (B-9b). Finalmente, Ao3 
se obtiene de (B-10c), quedando los coeficientes 


Ao3 = =09 > a 

Ar = pal 

Aj3 = pre (B-12) 
Aza = 209 

A>3 = E 


3 
d(x,y) = =pgy* — q (2 — LY + qpgay(z — 21). (B-13) 
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Comprobando el resultado. Veamos los resultados que se obtienen a partir de la función de Airy (B-13). Para las tensiones 
Tax, Se tiene 


04 04 3p9 
Oax = ya — P2Gs — PU9y = ya — P9U= =3 YE ay, (B-14) 


similar a lo que se obtendría de la ecuación (1-61), sustituyendo p = M/(2abL) y el momento de inercia de la sección plana 
rectangular por su valor /,. = 2ba? /3. Para las tensiones cortantes Ty Se obtiene 


09  3pg 


ZP PI E 
Toy 910y 207 Y ala —L), (B-15) 


similar a lo que se obtiene a partir de la ecuación (1-73) con h = 2a y el valor de T(x) dado por (1-62). Finalmente se 
obtiene, para las tensiones 0y, 


Po Po Pg, (2,2 
Oyy = qa 7 PTGw — PY9y = 53 — PYg= y ¿(a — y), (B-16) 
ecuación que no obtuvimos, pero que nos expresa la forma de las tensiones normales en la dirección vertical a que se hace 
referencia en la página 26. 
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